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ABSTRACT 


This report documents Phase I of a UNITECH research 
project directed toward development of a high speed automatic 
process to be used to match gun barrel signatures imparted to 
fired bullets. This project was contracted by NASA-JSC and de- 
fined in paragraph 3*2.1. a of the contract task description. 

An optical projection technique has been devised to produce 
and photograph a planar image of the entire signature, and the 
phototransparency produced is subjected to analysis using 
digital Fourier transform techniques. The success of this 
approach appears to be limited primarily by the accuracy of 
the photographic step since no significant processing limita- 
tions have been encountered. 
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1. INTRODUCTION 


In recent years, some of the most visible contribu- 
tions by the engineering community to the overall improvement 
in man's quality of life have been throuyh application of com- 
puter based automatic machinery to highly sophisticated and 
otherwise time-consuming tasks. Generally, these have boon 
tasks which had long been well understood, but for which con- 
temporary technology could not support execution of the tasks 
on the scale required. A good example of such a task is the 
Forensics application of ballistic slug identification. The 
generation of legal evidence by studying gun barrel signatures 
imparted to fired ballistic slugs, popularly termed "ballistics", 
is still considered to be a modern forensic procedure, but 
large-scale application of the technique has been prevented 
by the length and complexity of the process. Still basically 
a laboratory technique, the procedure is one which is per- 
formed manually by a trained technician and calls for moderate 
skill with a microscope and considerable personal judgment. 
Typically, ballistic analysis is not performed routinely 
because of limitations on time and trained personnel, in spite 
of the fact that most investigators of crimes involving fin 
arms eagerly seek ballistics assistance. Therefore, ballisti 
is a good example of a process which has not yet made the tran- 
sition from laboratory demonstration to production level appli- 
cation and can benefit significantly from high speed automation. 
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The central element in ballistics technology is the 
gun barrel signature. During manufacture, the bore of a gun 
barrel is drilled into a piece of solid steel and then finished, 
usually with rifling grooves added, by means of a series of 
reaming and polishing operations. Although in the more expen- 
sive weapons this polishing is very thorough, it is impossible 
to completely remove all surface roughness. When a bullet is 
fired through the barrel, the microscopic bumps and ridges 
along the bore literally scrape tiny grooves along the length 
of the bullet as it passes through, so that it collects a 
record of the bore topography. This record is incomplete, 
insomuch that it reveals little about the longitudinal position 
of bore surface features; instead it contains the composite 
of the circumferential or angular positions of bore features. 

The effect of rifling is to contribute to repeatability of 
the bore signature from one bullet to the next. The rifle 
grooves effectively lock the bullet in a fixed angular rela- 
tionship to the bore so that as the bullet moves along 
collecting markings, it cannot randomly turn about its axis 
which would destroy any record of angular relationships 
between features at opposite ends of the bore. In addition, 
the rifling process is one which leaves still more charac- 
teristic markings on the bore surface. In summary, it is the 
angular relationship of longitudinal marks on the bullet 
which constitutes the unique gun barrel signature and provides 
a basis for determining whether or not two bullets have been 
fired from the same gun barrel. 
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Ballistic analysis is currently performed optically. 
A weapon which is recovered during the course of an investi- 
gation is test fired into a capture box, usually fil. Led with 
cotton string or sometimes soap. The resulting test slugs are 
then assumed to * osfcss a representative signature of the gun 
barrel under consideration. One of the test slugs is mounted 
on one turret of a split image comparison microscope while an 
evi lence slug recovered from the scene of a crime is mounted 
on the other. Figure 1 is a photograph of the imago viewed 
by the ballistics laboratory technician. Each half of the 
split image field shows a magnified portion of the side of 
one of the slugs. The technician rotates one of the slugs 
with its turret until he sees a prominent feature. Ho then 
rotates the other slug in search of the same feature. When 
he finds it, he carefully aligns the prominent features in 
the two halves of the split image and then examines the 
positional alignment of the minor markings present on the two 
slugs. If most of the markings in his field of view appear 
to correspond on both slugs, the technician carefully rotates 
both slugs together, examining the entire 360 degrees of the 
signature. Inevitably there will be some marks which cannot 
be matched on both slugs, and it is this disparity which must 
be evaluated by the technician based upon previous experience. 
Should the evidence slug be damaged, as is often the case, 
the evaluation becomes still more difficult since the decision 
must be made from only a partial signature. 
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FIGURE 1 

PHOTOGRAPH OF GROOVE AREAS OF WO 
BULLETS DISPLAYED BY A COMPARISON MICROSCOPE 
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The object of the research effort documented heroin 
is not to develop an automatic system which can perform signa- 
ture comparisons ns accurntelv as humans can. Instead, it is 
hoped that a high speed system will \ developed which can 
accurately perform coarso comparisons to limit the number of 
bullets which must be examined by human technicians in order 
to provide a more widespread monitor of weapon related crimes. 
At present, ballistics can only handle the task of confirming 
or disputing the investigator's hypotheses. If a high capacity 
screening system were available, routine gun barrel signature 
checks could be performed on a city-, state-, or perhaps 
nationwide scale so that crimes related by a common weapon 
could bo detected automatically. This capability could elevate 
the role of ballistics from one of evidence support and confir- 
mation to one cf multiple offense detection. 

The next chapter, titled Problem Analysis, provides 
an overview of throe methods which have been proposed to com- 
pare signatures. One of these methods is chosen to bo the 
foundation iOr this research project, and the reasoning behind 
that selection is presented. Chapter Three is called Data 
Extraction and deals with those considerations affecting the 
spatial data sampling technique. Chapter Four describes the 
signal processing techniques applied to the signatures and 
presents the information obtained therefrom. The last chapter 
summarizes the progress made toward matching signatures and 
discusses limitations encountered and means by which these may 
be overcome. 
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The cjun barrel signature ns found on a bullet is a 
three-dimensional pattern, albeit a somewhat simplified one 
because it is restricted to the surface of a relatively simple 
geometric body. A further simplifying observation is that the 
pattern is situated on a single continuous surface. One of 
the limitations inherent in visual analysis of three-dimen- 
sional patterns is that only part of the pattern can be viewed 
at one time. It is this limitation that causes the balli3ti.es 
technician to turn the bullet to locate the salient features. 

At least from an optical standpoint, such patterns are much 
more easily recognized when they jro planar. Several tech- 
niques have been suggested for "unwrapping** the signature from 
the bullet. The most direct method is used by the profilo- 
meter, which consists of a stylus whose vertical travel is 
recorded as it follows the surface of a bullet which is turned 
beneath it. A lead screw moves the stylus along tne length of 
the bullet as the bullet turns so that the entire cylindrical 
surface can be probed. In many ways the profilometer resembles 
the original Edison phonograph. Another approach makes use of 
the same basic mechanical arrangement, but the bullet surface 
features are illuminated by a narrow, highly collimated laser 
beam, and the scattered light is detected by a photomultiplier. 
A third scheme produces a reflected optical image of the en- 
tire signature while avoiding the necessity for a rotating 
mechanism. Figure 2 is a cross-sectional diagram of a bullet 
centered in a conical reflector. This reflector projects an 
annuJar image of the bullet's cylindrical surface which can be 
recorded on film for subsequent digital processing. Figure 3 
is a photograph of a jacketed, hollcw-point bullet's image 
produced in this manner. 
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FIGURE 2 

DIAGRAM OF THE BULLET PHOTOGRAPHIC SYSTEM 


Tho third method is the one chosen for this project. 
It has the advantages that a planar image displaying tho entire 
signature is produced, and a photographic record is easily ob- 
tained which permits both digital and optical processing. 
Digital processing begins at this point by scanning a photo- 
transparency of tho signature imago with a microdensitometer. 
The microdensitomctcr measures optical density of the photo- 
transparency points defined by a two-dimensional X-Y grid. A 
density value is generated and stored on digital magnetic tape 
for each X-Y coordinate pair. Signatures digitized in this 
manner can then bo compared by cross-correlation or matched 
filter techniques. 
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Tho phototranaparcncy can also bo u^ud to match siq- 
natures in a totally optical system. This approach exploits 
tho capability of lenses to produce two-dimensional Fourier 
transformations. Figure 4 is a diagram depicting an optical 

system which will produce tho two-dimensional Fourier trans- 

* 

form of a phototransparency. Goodman r l] explains that a 
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FIGURE 4 

DIAGRAM OF AN OPTICAL FOURIER PROCESSING SYSTEM 


converging (convex) lens produces the exact two-dimensional 
Fourier transform, in its rear foca 1 plane, of a two-dimensional 
object image located in its front focal plane. The object imago 
is assumed to be illuminated by a monochromatic plane wave. In 
reality, the range of the transform is restricted to the primary 
lobe of the Fresnel diffraction pattern for the lens aperture, 
but this limitation can bo easily accommodated if the lens can 
be made large when compared with the object image. If a diver- 
ging coherent reference beam is directed onto a photographic 


♦Bracketed numbers correspond to reference numbers in the 
bibliography. 
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emulsion located in the transform plana in Figure 4, a holo- 
graphic matched filter for the object transparency will bo 
recorded. For optical matched filter processing, the filter 
phototransparency is replaced in the transform plane of the 
same opcical system but no reference beam is used. If the ob- 
ject transparency contains the image to wnich the holographic 
filter is matched then a reference beam will be holographically 
reconstructed which is focused on the focal point of the ori- 
ginal diverging reference beam. A photo detector located at 
that point serves to signal that a match has occurred. 

Still another processing method is available for 
transparencies which is an optical-digital hybrid. In this 
method, the converging lens is used to generate the Fourier 
transform just as in the totally optical system, but it projects 
the spectrum onto a photodiode detector array situated in the 
roar focal plane. Each diode in the array produces a voltage 
proportional to the intensity of light striking it, so that if 
the array outputs are sampled and digitized, a digital two- 
dimensional power spectrum is obtained. While a power spectrum 
cannot be used to perform cross-correlation, it nught be useful 
to perform coarse signature screening of the basis of power 
spectral comparison alone. If a digital data library of stored 
normalized power spectra were available in a digital computer 
interfaced to the diode array, then it would be a simple matter 
to insert a transparency of a new signature into the optical 
i ssembly and command the computer to search its library for a 
matching spectrum. 

From the foregoing presentation of alternatives, it 
would seem that the optical systems should be preferred from 
the standjx^int of simplicity. As always, however, there are 
problems inherent with optical systems which theory does not 
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address. The primary areas of concern in this application can 
all bo grouped into the category of mechanical alignment. 

There is degradation in optical systems in yenoral trom imper- 
fections in lens and focal plane gnome? tries. Thio specific 
application requires in addition that careful attention be paid 
to image positioning and size on the transparencies. A perfect 
match of signatures requires not only that the signature charac- 
teristics be matched and clearly di scernable, but al no that the 
signature image on the object transparency bo identical in size, 
position, and orientation to the image. This latter require- 
ment means that in the case of the purely optical processor, 
the object transparency must bo centered on the lens axis and 
rotated through 360 degrees to determine the best signature 
correspondence orientation. All of the above i msiderations 
place accuracy requirements upon either the photographic step 
or the optical system operator's placement of the transparency, 
thereby defeating our intention to develop a fast efficient 
system. 


The all-digital system is somewhat cumbersome when 
the microdcnsitometor is remembered, but it provides the most 
accurate and flexible processing available. Since it is not 
clear initially what degree of resolution will be required to 
match signatures with the accuracy desired, it is prudent to 
perform the initial research with as much accuracy cind flexi- 
bility as possible. Therefore, the research described heroin 
is done with the digital system with intent to determine per- 
formance capability of an "optimum" system. This optimum or 
baseline system will be used in future work to perform statis- 
tical measurements as functions of several parameters. From 
those measurements, it should be possible to establish minimum 
performance requirements for use in determining feasibility of 
the optical techniques. 
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III. DATA KM TRACT in;; 


Digital processing of tho ballistics signatures begins 
by sampling and digitizing tho phototrnnspnrencios made using 
the conical reflector. These transparencies arc "negatives", 
one of which was used to produce the positive print in Figure 3. 
Since the signature information is found in a two-dimensional 
pattern, two-dimensional sampling is required. The scanning 
microdonci tome ter used to sample and digitize the signature 
transparencies ii* programmed to perform two-dimensional sampling. 
It consists of a glass table which can be positioned in X and Y 
coordinates with a resolution of 1 x 10 ’ meter (1 uM) and an 
adjustable aperture collimated light source which projects a 
beam through the table in’O a detector. In operation, the 
transparency is measured by the detector at every sample point. 
The detector output is digitized by an analog-to-digital (A to D) 
converter and the resultant digital density values for each 
point are stored on magnetic tape together with tho coordinate 
values . 


The location of the planar sample function on tho 
transparency is established by defining the origin for the 
microdensi tome ter ' s Cartesian coordinate system such that the 
entire bullet image is just contained in the first quadrant, 
as limited by the maximum coordinate values. Tho maximum 
values are the same for both X and Y coordinates and are de- 
fined primarily by tho sample spacing. The maximum allowable 
spacing must be determined from the data bandwidth to reduce 
data storage requirements . Maximum sample spacing can be 
determined graphically in this application with the aid of a 
calibrated microscope and a representative group of trans- 
parencies. The signature negatives are surveyed to establish 
tho minimum width of a signature line. This measurement must 
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ho made at tho inner. -most part of the signature. The reason is 
clear when the image mapping mechanism is considered. Figures 
2 and 3 demonstrate that the gun barrel signature as embodied 
in tho cylindrical sur ce of the bullet is mapped into an 
annular imago by the truncated conical reflector. If tho bullet 
is positioned point outward as in Figure 2, the base will be 
mapped along the inner edge of tho annulus while the point 
appears as the vaguely defined outer edge. Signature lines of 
constant width along the length of the bullet are mapped into 
radial linos of constant angular width in tho annular image. 
Therefore, a signature line in the image is narrowest at the 
end closest to tho center of the reflector, which is the inner 
edge of the annulus and corresponds to the base of the bullet. 

A sample spacing chosen to capture tho smallest signature lines 
in this area will certainly suffice at all points for all lines 
in the signature image. The criterion for choosing a maximum 
sample spacing can bo understood with the aid of Figure 5. 

Figure 5a is a simplified diagram of the sample locus foi a 



FIGURE 5 

IMAGE SAMPLING DIAGRAM 


single scan along cho inner edge of the image annulus displayed 
in an area wiicre it is normal to two dark minimum width signa- 
ture lines separated by a bright line. The objective is to 
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space sample points, indicated by dots on tho locus, so that 
none of those linos are skipped; and this is accomplished as 
indicated for tho worst rase if the points are spaced slightly 
closer than the width of tho narrowest lino. This graphic 
criterion is equivalent to the Nyquist criterion, or what 
Braccwell [2] calls "critical sampling" as con bo seen in 
Figure 5b. Tho sinusoid vn part b represents the transmitted 
light intensity, or Z amplitude, of the lines in part a with 
the corresponding sample points indicated. The Nyquist crite- 
rion, simply stated, says that the sampling frequency must b'* 
sligntly greater than twice tho maximum frequency present in 
the analog signal, and clearly it is met. 

The foregoing analysis servos to ascertain that the 
desired signature bandwidth will be preserved in the sample 
data, but it does not assure that the imago on the transparency 
is \ mdlimitod below the Nyquist frequency'', as it must bo to 
prevent aliasing of high frequency components into the signa- 
ture band. The microdensitometer can perform a low-pass fil- 
tering function on the image data before it is sampled, however, 
by averaging with the adjustable beam aperture. If a square 
aperture is selected whose width is equal to the sampling inter- 
val, then each sample will represent the average of all imago 
points in the beam and these averaging areas will be exactly 
contiguous. This scheme efficiently serves as an anti- 
aliasing filter. 

After the bullet image is X-Y sampled and digitized 
it will be used in several ways for comparison with other 
images, notably by correlation. If the two-dimensional corre- 
lation process is pictured graphically as sliding one image 
around over another while plotting the sum of products of all 
pairs of sample points for each relative position of the images. 
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it is obvious that some serious problems exist. The most ob- 
vious is that there are too many mul tipi icat ions required. A 
typical sample array contains 1024 x 1024 words, an over- 
whelming size for correlation processing. A still more serious 
problem appears when oi.o considers that two-dimensional cross* 
correlograms are functions of X and Y lags only, containing no 
information about relative image rotational alignment. Although 
it is practical to center the images in X and Y before they are 
sampled, l hero remains the problem of performing two-dimensional 
correlation as a function of rotational lag on X-Y samples 
arrays. Such a process would require interpolation to locate 
corresponding samples in the two arrays as one is rotated. 

If the sample spacing criterion is recalled, the 
bullet image is oversampled everywhere except along the inner 
edge of the annulus, which indicates that possibly the data 
volume could bo reduced. Figure Ga is a diagram of the bullet 
located in a cylindrical coordinate system, and part b dia- 
grammatically shows the bullet with cylindrical coordinates as 
mapped by the conical mirror. It is apparent from Figure Gb 
that the digitized X-Y image can be resampled in h-0 coordi- 
nates so that a sample array resembling Figure 6c results 
which is in essence the signature unwrapped from the bullet. 

The h and 0 sample spacing is chosen equal to the original X-Y 
increments at the inner edge of the annulus where is crosses 
the X-tixis so that the maximum sample spacing is obtained through- 
out the image. The image center is determined manually on the 
microdensitometer arid routinely defined to be 512, 512 in the 
X-Y coordinate system. Resampling is accompli shed by defining 
a "ring" of 0-points for each incremental h value and inter- 
polating the sample values at those points from the X-Y sample 
array. Since only the side of the bullet contains signature 
information, the h-axis can be truncated to eliminate the crimp 
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ring around the bullet near tho point as well as the joint 
itself. This is a very useful capability since it is common to 
recover bullets damaged on the front end. Figure 7 presents 
a useful illustration of tho resampling process. It is a con- 
tour plot of a small portion of a bullet image interpolated 
from the X-Y sampled data as digitized by the microdonsi tome ter . 
The contour linen enclose local amplitude peaks, or bright areas 
from the planar bullet image. Superimposed on tho contour plot 
is a constant-h resampling locus corresponding to a ring around 
the bullet approximately halfway between the base and the crimp 
band. Figure 0 is a plot of the amplitude samples interpolated 
along the same locus shown in Figure 7 and to the same scale. 
These now samples are representative of tho data present in tho 
rings after resampling. 

Another consideration encountered when extracting data 
with a microdonsitometor from images recorded on photographic 
film is the transfer function of tho film. Kodak [3*1 explains 
that within the linear jiortion cf the film characteristics curve 
photographic emulsions darken in jiroportion to the logarithm^ 
cf the intensity of incident light. Tho exact relationship is 
expressed as 


D “ 

whore P is optical density and I is intensity. A microdonsi - 
tometer of the type used in this research commonly reads out 
this digital density as it scans the transparency. Since den- 
sity represents the logarithm of the original image intensity, 
considerable dynamic range compression is present. This com- 
pression reduces accuracy in all of the contemplated processing 
techniques because of its non-linear filtering effect. To avoid 
non-1 i neari ty density is converted back to intensity, but the 
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FIGURE 7 

CONTOUR PLOT OF A PARTIAL X-Y DATA ARRAY 
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conversion must bo per formed before resampling since linear 
interpolation of density values does not produce linearly inter- 
polated intensities. 

Another observation can be made contrasting the X-Y 
and h-9 sampled data representations. The X-Y format of the 
imago data requires two-dimensional processing for all signa- 
ture comparison schemes because of the annul ur shape of the 
signal tire region in the X-Y plane. In the h-0 plane, however, 
the signature occupies a rectangular region. Moreover, reflec- 
ting on Figure 6c, it seems likely that sufficient information 
to compare signatures can be found in the relative positions of 
lines of the signature, described as a function of 0 alone. 
Resampling conceivably reduces the two-dimensional problem to 
one dimension. 

In contemplating the one-dimensional approach, addi- 
tional raw data preparation steps appear advisable. The contour 
plot in Figure 7 demonstrates that signature lines are not of 
constant amplitude along their lengths, and in fact, some even 
approacli zero amplitude in places. This problem suggests that 
heavy averaging in the h dimension might be useful. This would 
remove h-depondent high frequency information, but the onc- 
iimcnsional scheme presumes that very little h-dependont signa- 
ture information exists. If very many rings are to be averaged 
together, however, some deskewing must bo done to compensate 
for the twist imparted by rifling to the signature. This 
rifling twist prevents the signature lines from being strictly 
parallel to the h-axis, but the skew is linear with h and easily 
measured. Combined with the other gross characteristics such 
as caliber and number of grooves skew or twist makes an excel- 
lent signature library categorization clue. 
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Admittedly, ona-dimonsional analyois does not provide 
intermediate outputs that aro directly comparable to those from 
the optica] processors being considered. Of particular inte rest 
would be the fowor spectrum output for comparison with the op* 
tical power spoctrum processor output, but a basis for perfor- 
mance prediction will still exist. The one-diinonaional power 
spectrum should provide a worst-case test for power spectrum 
based processors, and should it. prove to be a useful signature 
discriminator, the optical version should work at loast ar. well. 



IV . KXPERIHgHM KTS 


Tho research conducted for this project is part of a 
joint effort f undcd by tno National Aeronautics and Space 
Administration (NASA). Tho photographic signature recording 
was done by Grumman Aerospace Corporation using t.ho roflectivo 
con<* assembly previously described. Only tho digital processing 
approach was to bo invootiga ted by UNITECH in Phase 1 of tho 
project, using a number of phoLotransparencies digitized with 
12-bit resolution on Perkin-Elmor model 1010A microdensi tome tor 
systems located both at NASA's Johnson Space Center (JSC) and 
at the University of Texas at Austin. Signature analysis v/as 
performed with a microscope equipped with a calibrated ocular 
micromotor on tho 35 mm transparencies# from which it was deter- 
mined that a square sample grid of 1024 x 1024 points would be 
required to assure adequate imago resolution. After sampling, 
tho digitised image denwition wore converted to linear intensity 
values and resampled in h and 0 increments chosen to maintain 
equivalent image resolution. Tae resulting imago sample array 
contains 58 rings spaced at uniform h increments up to but 
excluding the bullet crimp band. Each ring contains 1450 sam- 
ples uniformly spaced in 0. Tho exact number of 6 increments 
chosen was determined from processing sj>ecd considerations in- 
herent in the transform processing next to be introduced. 

One of the early concerns about tho digitized image 
was the bandwidth of the signature. A quick check can be made 
using the power spectrum of one or more rings since little or 
no high frequency signature information is expected to exist 
as a function of h. The paw».!r spectrum can bo obtained from 
the Fourier trans form of tho ring data function, since 
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r ( f ; « k x h * 


whore 


R *■ Ffr] , 

* 

and R is tho complex conjugate of R. This metnod for obtaining 
tho power spectrum is practical because of tho speed of current 
Fast Fourier Transform (FIT) algorithms. Tho fastest ones 
available require that the number of input sample's bo an inte- 
ger power of two, but in this application r.uch a const raint was 
judged hazardous, at least until tho useful bandwidth of tho 
signature could be determined. The minimum radius of the sig- 
nature annulus is sampled by the microdonsi tomotcr is 231 incre- 
ments long which means that the first ring has a circumference 
of 1451 increments. Rather than undcrsample at 1024 points per 
ring, it was decided to use a mixed-radix ITT and tho one des- 
cribed by Singleton [4] war. selected. Its speed depends upon 
tho prime .‘.'actors of the array length integer, and 14 38 is tho 
closest integer to 1451 permitting efficient FFT operation, thus 
1458 became the ring array size. Figure 9 is an unnormalized 
power spectrum generated from a representative sample ring 
after transformation using the FFT. ic was plotted as a single- 
sided spectrum and extends from a spatial frequency of 1 to 729 
cycles per revolution (CTR) which is the folding or Nyquist 
frequency. The dc component was removed, but no ring power 
spectrum averaging or other smoothing was done since tho only 
objective was to roughly determino tho signature bandwidth. 
Obviously .it: falls well within the 729 CPR band, and it seems 
doubtful that any important components exist above 304 CPR or 
half of tho current bandwidth. 


Power (Linear scale) 



FIGURE 9 

UNNORMALIZED POWER SPECTRUM PLOT 
OF A SINGLE :.'ATA RING 
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Given the above information it is reasonable to re- 
duce the ring array si ze to speed up processing. Since the 
original bandwidth war. twice that deemed necessary, it was 
decided to a* erago and decimat e the ring data camples by two. 
Averaging pairs of samples effectively halves the bandwidth to 
prevent aliasing so that only half of the original number of 
ring data samples must be used, specifically the averages of the 
original pairs. 

The next technique to bo investigated was the previ- 
ously mentioned averaging of data rings as an aid to one-dimen- 
sional processing, using 9-dcp.ndent information alone. Prior 
to averaging a large number of rings, an effort was made to 
measure the skew present as a function of h. The method used 
consists of averaging the first and last five rings to provide 
two separate rings representative of signature data at opposite 
ends of the bullet, and then cross-correlating these rings to 
measure the rotational lag or skew between fc v *em. Again making 
use of the FFT, the cross-correl ogram is expressed as the in- 
verse transform of the complex cross -power spectrum of the two 
rings . 


F 1 [R^ x R^, 1 « cr 12 ( t) 

Using Leo's notation [51, a positive lag, t# corresponds to a 
positive rotation of ring 2 with respect to ring 1. Figure 1C 
is the corrclogriim which was generated for bullet A-4 except 
that the DC component was removed from the cross-power function. 
Two particularly disturbing bits of information appear in this 
plot. First, there *s very little skew measurable between the 
rings as indicated by the location of the peak near the origin 
on the t axis. Second, the next highest peak is comparable in 
size to the main peak. Those observations prompted a closer 
examination of the ring cample function to find the explanation. 


Cross Power (Linear scale' 
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FIGURE 10 

UNNORMALIZED CORRELOGRAM OF WIDELY-SEPARATED 
FIVE-RING AVERAGES FROM BULLET A-4 
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Figure 11 is an intensity plot of the lower of the 
two five-ring averages. The largo amplitude spikes correspond 
to the brightest signature lines in the image which emanate from 
the edges of the rifle grooves. Since there are five grooves, 
there are ten edges which should produce ten equally spaced 
spikes. Several of them are missing, and a re-examination of 
Figure 3 shows them to be missing in the image, probably due to 
non-uniform lighting when the transparency was made. The ex- 
treme amplitude disparity between the groove boundary spikes and 
the remaining signature structure suggests thol the correlogram 
in Figure 10 is dominated by spikes. The effect, of this domi- 
nation is that if the correlation peak amplitude were to be used 
as an indication of signature similarity, then any two signa- 
tures with the same rifle groove pattern would be judged vir- 
tually identical. This conclusion is supported by the obser- 
vation in Figure 10 that a very high correlation peak occurs at 
a lag vaJue equivalent to a misalignment of onu groove width. 
Fortunately, at least one simple solution exists for this pro- 
blem. If the spikes are limited or clipped at a very low ampli- 
tude, the secondary peaks should disappear from the correlogram. 
This solution is supported by Figure 12 which displays the 
corrclogx-nm generated using clipped ring data. 

The cause of the poor skew indication in Figures 10 
and 12 was not clear at this point in the research. The ampli- 
tude of the norma] ized peak is loss than 0.5 instead of almost 
1.0 as had been do ;ired but some correlation reduction might bc- 
expcctod if the signature structure is not totally independent 
of h as this one-dimensional analysis assumes. 

The next stop performed was cross-correlation of 
averaged rings from signatures of two different bullets fired 
from the same gun. The normalized peak value of this correlo- 
gram was expected to give some indication of the degree to which 
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FIGURE 11 


INTENSITY PLOT OF A FIVE-RING AVERAGE 
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tost fired bullet signatures would match and thus suggest where 
the match detection threshold might bo set. Because the skew 
measurement had boon questioned# ensembles of the individual 
rings to bo averaged were plotted for each signature image. 

Tho ring plots wore stacked one above the other to determine 
how much signature structure might bo lost in tho ten-ring 
average, and whether or not uny deskewing would be required. 

Those plots arc found in Figure 13. In these plots several ob- 
servations are important. First# although good comparison of 
the tv/o signature ensembles is not expected at zero lag, it 
appears that they do not compare at any lag value. It had been 
expected that at least the unclipped groove spikes would match 
at ten different relative orientations# but they do not. When 
measurements were made of tho spacing between rpikes in each of 
the bullet signatures# and allowing lor missing spikes# they 
were found to be vary non-uniform. Furthermore# close inspec- 
tion of the skew across tho ensemble for each spike shows that 
the skew is non-uniform among the groove spikes within the same 
signature, which probably explains tho poor performance of the 
skew corrologram. This distortion of gross image structural 
features can always be expected to ruin correlation processing, 
and indeed in Figure 14 the cross-eorrelogrnm using these 
averaged, but clipped ring ensembles has no primary peak at all. 

A close examination of photo enlargements of tho 
transparencies from bullets A2 and 7*4 showed that the angular 
width of rifle grooves in tho images is irregular, and moreover, 
that tho rifle groove pattern on one photo cannot be matched 
with the pattern on Die other because of some angular distortion . 
Furthermore, those are signature patterns from two bullets which 
had been found to match closely when visually observed under a 
comparison microscope Based on these observations, it was 
postulated that some misalignment in the photographic stop 
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(a) 

TEN-RING ENSEMBLE PLOT OF INTENSITY DATA FROM BULLET A-4 
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(b) 

TEN-RING ENSEMBLE PLOT OF INTENSITY DATA FROM BULLET A- 2 


FIGURE 13 
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superimposes an angular distortion function onto tho signature 
image. Tho orientation of this function relative to tho ro- 
tational position of tho bullet varies from one bullet to tho 
next since no attempt is made to rotationally j>osition each 
bullet in tho cone tho came way. If this explanation wore in- 
deed the cause of the problem, it was reasoned that two photo- 
graphs made of tho same bullet without moving anything should 
exhibit signature images distorted in exactly the same manner, 
and cross-correlating tho two should produce a well defined peak 
at zero lag. Fortunately such a pair of photo transparencies had 
been made as a test of exposure settings, and the resulting 
croas-correlogram is presented in Figure 15. Tho fact that the 
expected peak is present vindicates all of the digital pro- 
cessing steps, includin the re sampling process which would have 
boon the next logical suspect. However, this discovery halted 
further signal processing analysis pending resolution of tho 
photographic problems. Subsequent visual observation of a 
bullet in the conical reflector has demonstrated that the type 
of distortion encountered in the transparencies is produced when 
the camera lena ia not centered on the . < , 

no improved transparencies have been made at this time. 
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FIGURE 15 

NORMALIZED CROSS-CORRELOGRAM OF TEN-RING AVERAGES 
FROM TRANSPARENCIES 7 A AND 9A 



V . SUMM A RY AND CONCLUS JONS 


Research reported heroin has applied some common digi- 
tal signal processing techniques to the problem of recognition 
ol gun barrel signatures imparted to fired bullets. Primary 
tools utilized have been spectral analysis and correlation 
analysis, both of which were implemented with the aid of tho 
Fast Fourier Transform (FFT) . Signature extraction was accom- 
plished, with as yet only moderate success, by photographing 
the signature image as projected onto a plane by a conical re- 
flector. Since both digital and optical processing can bo 
applied to the resulting photo transparencies, it was hoped that 
tho highly flexible digital approach would provide a reference 
level of performance which might be used to evaluate simpler or 
fnslor optical or hybrid systems. Unfortunately only a cursory 
ovalui n ol the digital techniques mentioned has been possible 
due to alterations in tho optical components of the photographic 
system which severely distorted the signature images* 

Tho processing techniques described heroin have demon- 
strated encoui aging capability. With little or no frequency 
domain filtering it is possible to produce corrolograms with 
well defined peaks. Tho power spectra of representative signa- 
tures are well behaved, and the sampling scheme developed per- 
mits ensemble averaging of power spectra to support meaningful 
spectral analysis. As soon as a practical method is developed 
to produce accurate signature transparencies, this effort should 
prove successful. 
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